On the surface of glasses 
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Dynamics near the surface of glasses is generally much faster than in the bulk. Neglecting static 
perturbations of structure at the surface, we use random first order transition theory to show the 
free energy barrier for activated motion near a free surface should be half that of the bulk at the 
same temperature. The increased mobility allows the surface layers to descend much further on the 
energy landscape than the bulk ordinarily does. The simplified RFOT calculation however predicts 
a limiting value for the configurational entropy a vapor deposited glass may reach as a function of 
deposition rate. We sketch how mode coupling effects extend the excess free surface mobility into 
the bulk so that the glass transition temperature is measurably perturbed at depths greater than 
the naive length scale of dynamic cooperativity. 



Paying attention to the properties of interfaces has, 
since the time of Van der Waals, helped to clarify our 
understanding of bulk phase transitions [T] . Studying the 
mobility at the surface of glasses and supercooled liq- 
uids has the potential to be equally enlightening about 
the glass transition [2 3 . Theories of the glass transition 
predict a growing length scale of correlated dynamics as 
the glass transition is approached [H IID E| • This growth, 
according to the most successful theories, is quite mod- 
est. It has only in recent years been widely acknowledged 
that there is such a slowly growing length scale in bulk 
glasses Ej. The dynamics of a glass or supercooled liq- 
uid should be perturbed within a few correlation lengths 
of its surface. Many experiments do show significant per- 
turbations of the glassy dynamics at free surfaces or in 
confined spaces [TO] ■ The picture that emerges from these 
experiments is, at present, still somewhat confused. In 
general, mobility seems to be increased at a free surface, 
although sometimes a diminished mobility has been ob- 
served. At interfaces with solids both increases and de- 
creases in mobility have been reported |lll 112] . It seems 
likely that this complex set of behaviors reflects the fact 
that the dynamics of the liquid can be strongly influ- 
enced by the liquid's static structure which also will be 
perturbed by the interface — in the extremes, partial 
crystallization can occur at a free interface while a dry- 
ing layer may sometimes insulate a confined fluid from its 
solid surroundings. While noting these complications, we 
feel it is nevertheless worthwhile to present a simplified 
view of glass and supercooled liquid surfaces in the con- 
text of the random first order transition theory of the 
glass transition. 

The random first order transition (RFOT) theory is 
a constructive approach to structural glass dynamics 
that has explained quantitatively numerous bulk glass 
phenomena^. It was recognized very early on that grow- 
ing length scales should be associated with an ideal glass 
transition as envisaged in RFOT theory [T3| H4], Indeed 
some of the early experiments on confined supercooled 



liquids carried out by Jonas [15] were motivated by a de- 
sire to test these expectations. It turns out present day 
RFOT theory for molecular liquids can make somewhat 
more definite predictions than was done in those early 
days, at least for an idealized interface which can be taken 
to have no change in static structure from the bulk. We 
show for this idealization that the maximum mobility at 
a completely free interface is related in a very simple 
way to the bulk mobility. In terms of relaxation times, 
the RFOT result for the surface relaxation time is simply 
Tsurf = ^/TQTbuik- We also show that RFOT theory sug- 
gests that as measured by an effective local glass transi- 
tion temperature T l ° cal the influence of the interface can 
appear to be rather far-reaching into the bulk, consistent 
with some experiments. 

At the smallest length scales the dynamic arrest of 
glasses results from the cage effect in which a particle's 
motion is constrained by the presence of its neighbors. 
In mean field theory this leads to a friction crisis at the 
mode coupling critical temperature and the emergence 
of an ensemble of aperiodic crystal structures having an 
extensive configurational entropy s c per particle. Below 
the dynamical transition at T c any large scale motion 
that takes place is necessarily collaborative. Near a free 
surface this picture is modified. Since surface particles 
feel a weaker structural cage, essentially only on the in- 
ner side, they would go through a dynamic arrest at a 
lower temperature and remain more mobile below the 
bulk glass transition temperature. The collaborative dy- 
namics propagates the enhanced mobility at the surface 
some distance into the bulk. This depth would be deter- 
mined by the length scale of cooperativity. 

In the mean field limit (made precise in Kac models by 
Franz [TB]) two cooperative length scales can be defined 
for random first order transition theory [TJ]. One scale, 
£mct, is directly related to the dynamical transition at 
T c which resembles a spinodal[7], the other is related to 
the size of regions that rearrange by activated motions 
and scales with the configurational entropy s c [l7 . The 
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latter scale diverges at the Kauzmann temperature with 
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an exponent ^rfot ~ cls c ■ Within the sharp, thin 
wall approximation the RFOT theory yields a numerical 
coefficient a that appears consistent with experimental 
values for molecular glasses 17 \. The scaling exponent 
used to fit inferred lengths from \i is a lso n °t inconsis- 
tent with RFOT theory, although a larger value provides 
a better fit [SI E] ■ For molecular liquids governed by short 
range interactions the dynamical (mode coupling) length 
would not actually diverge at T c because of a dynami- 
cal cutoff from the activated events and will not be too 
different from £,rfot- 

We first review the theory of bulk activated 
dynamics[l]. Below the mode coupling theory (MCT) 
transition, dynamics takes place on a rugged free en- 
ergy landscape dHHSl US]. Particle motion occurs locally 
through transitional hops between metastable structural 
states which resemble, in many ways, nucleation pro- 
cesses between different aperiodic crystal structures [HI 
117] . The magnitude of the free energy barrier can be 
found by a competition between the entropic cost of re- 
maining confined to one minimum free energy structure 
and a mismatch free energy penalty, a, for having two 
mean field solutions adjacent to each other but in dis- 
tinct structural states. The free energy profile governing 
nucleation of a spherical cooperative region of radius r is 

F(r) = 4 7 rr 2 (r /r) 1 / 2 <7Q - ^TTr 3 n Ts c (1) 

The surface penalty a{r) = (ro/r) 1 / 2 ^ scales with r 
due to renormalization effects of wetting by the diverse 
set of structures [T7| [5D]. The r -1 / 2 scaling holds suffi- 
ciently close to the ideal glass transition. At short length 
scales (Jo can be obtained from a crude density functional 
calculation giving a = ^fisTr^ 2 In where <Il ~ 0.1 

particle spacings is the Lindemann length, no = r^ 3 is 
the density of the liquid. The configurational entropy, s c , 
measures the number of available structural states. The 
resulting free energy barrier for reconfiguration events in 
bulk is 
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FIG. 1: In the bulk, at low temperatures, activated motion 
occurs within roughly spherical regions. Near a free surface 
a rearranging hemispherical region feels no mismatch penalty 
on interface along the free surface leading to much faster dy- 
namics than the bulk. 

This length scale increases with decreasing temperature, 
reaching a universal value of about 5 inter-particle spac- 
ings at the glass transition temperature corresponding to 
a one hour relaxation time. 

The RFOT analysis of bulk activated dynamics can be 
easily modified (see figure [T]) to treat motions near a free 
surface. Near a completely free surface with no struc- 
tural modifications, the transition state would rearrange 
a region of hemispherical shape in order to minimize the 
surface area subject to the mismatch penalty. This shape 
can be shown to be appropriate when wetting effects are 
neglected [31]. The flat face of the hemisphere lies along 
the free surface and is assumed not to contribute at all 
to the mismatch penalty. The free energy profile for the 
hemispherical rearranging region (including wetting) be- 
comes 



F surf (r) = 2irr 2 (r /r) 1 / 2 l Jo ~ -irr 3 n Ts c . (4) 

Since both terms in the expression are halved from the 
equation for bulk dynamics we see that the size corre- 
sponding to the transition state for activated rearrange- 
ment remains unchanged from that for the bulk, but the 
resulting free energy barrier is reduced by a factor of 2, 
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The relaxation time r a = T e Ft / kBT diverges at the 
Kauzmann temperature (s c (Tk) — > 0) and thus follows 
the Vogel-Fulcher behavior in the deeply supercooled re- 
gion r a ~ e B /( T - T K) _ The length scale of an activated 
event also follows from the free energy profile 
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This seemingly innocuous change leads to dramatically 
faster relaxation at free surfaces than occurs in the bulk. 
Molecular motion near the surface of a glass would still 
be detectable even if the bulk were essentially frozen. 
We see that very near the surface, the relaxation time is 
related very simply to the bulk value, 
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According to RFOT, on laboratory time scales, the 
ideal surface layer will be able to "hike" (using the color- 
ful expression of Ediger [2l]) much further down the free 
energy landscape than expected, reaching lower energy, 
more stable configurations than the bulk is likely to find. 
The surface motion will freeze out, i.e. go through its 
own glass transition (Tg Ur f) only when r sur f becomes 
larger than the laboratory time scale. For the same lab- 
oratory time scale, according to equation [5] this will hap- 
pen when the configurational entropy is half the char- 
acteristic configurational entropy of the bulk glass tran- 
sition. The bulk glass transition entropy is very nearly 
material independent [22j [23], and likewise the free sur- 
face glass transition will occur at a characteristic bulk 
entropy value. For a one hour time scale this character- 
istic entropy is 



yielded 9k ~ 0.4 for both liquids studied |21_. It would 
take 10 4 years to to make these highly stable glasses via 
traditional bulk cooling, while 10 9 years, roughly the age 
of the universe, would be required for a glass with the 
limiting stability of 9 k = 1/2. Unless the surface struc- 
ture is greatly modified leading to a locally smaller Tk, 
vapor deposition should not yield glasses with stability 
much greater than the naive RFOT limiting value. 

The linear extrapolation for s c (T) is valid near Tk, 
but for some liquids, at higher temperatures a better ap- 
proximation is |24j s c = Soo(l — Tk/T). Recognizing this 
feature does add some material dependence to the limit- 
ing stability 
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Ediger and co-workers |21j have recently shown how to 
construct a microscopic glass sample via vapor deposi- 
tion reaching much lower in the energy landscape than 
was possible by the usual bulk cooling of a liquid. They 
argue this ultra-stable glass arises from the free surface's 
excess mobility. For slow deposition on a cold substrate 
the growing glass, owing to the higher mobility, can more 
thoroughly rearrange near the free surface to form ex- 
ceedingly stable structures. They assign a Active tem- 
perature to their sample, Tf, defined as the temperature 
where a glass of similar stability, but created by bulk 
cooling would fall out of equilibrium. Tf can be used to 
measure the sample's position on the free energy land- 
scape relative to an ordinary glass through the parame- 
ter, 



Fragile liquids still have a limit near 1/2, but strong liq- 
uids have the possibility of descending further in terms 
of temperature (but not in terms of configurational en- 
tropy) than fragile liquids. The fragilities for the two 
glasses so far studied by Ediger are similar and would 
yield 9 K < 0.57(IMC) and 9 K < 0.58(TNB) for the 
nanometer per hour deposition rate. 

A key prediction of this analysis is that the Active tem- 
perature is a logarithmic function of the deposition rate. 
More precisely the configurational entropy at the Active 
temperature is a universal function of the logarithmic 
deposition rate. The limiting Active temperature given 
above assumes a deposition rate of one correlation length 
per hour. More generally, for a deposition rate k the mo- 
bile surface layer of depth £ will be equilibrated on a 
time scale which is related to the Active temper- 
ature through the relationship £/fc = TQe F ( T /)/ fcfiT / = 

Toe A'/s c {T f )_ 
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This parameter increases with the stability of the glass, 
reaching one for an ideal glass. In principle if there were 
in fact no entropy crisis, 9k could exceed one. 9k can 
be parametrized by the conAguration entropy, an equiv- 
alent measure of structural stability. Using the linear 
extrapolation vanishing at Tk valid for low temperature 
(T < T g ) the RFOT free surface mobility would yield, 
for a rate of deposition equal to one correlation length 
per hour 
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As the enhanced surface mobility freezes out at Tg Ur \ 
this temperature should be seen as a limiting value for 
Tf leading to the inequality above. Ediger's experiments 
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We compared to Ediger's experimental results using 
the linear representation for the conAgurational entropy, 
s c (T» = AC p (T g )(Tf- T K )/T K . Along with A' ob- 
tained from equation pi we take £ rs lnm, r ~ 10 _12 s. 
ACp(Tg) ~ 2.6fc_B is the heat capacity jump at the glass 
transition for IMC in Boltzmann units per bead [23 . The 
theoretical results along with the measured values ob- 
tained by Ediger et al.[2T] are plotted in Agure [2] Our 
results are consistent with the experimental data for the 
deposition rates tested. The theory correctly predicts the 
trend for slower deposition rates. Despite the good agree- 
ment this is a zeroth order approximation in which we 
have assumed equilibrium relaxation, strictly true only 
for T su b s t ra t e — Tf . 

It would seem natural to say that the enhanced mobil- 
ity at the surface would penetrate into the bulk at least 
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We have shown that A sur f can be as low as -^Abuik, 



and 
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FIG. 2: Fictive temperatures vs deposition rate for the glass 
former IMC. Data for the deposition experiment was taken 
from reference [5T]. The same data are shown in the inset 
plotted as s c (Tf) in Boltzmann units per bead vs deposition 
rate. 



on the length scales of the bulk dynamical heterogene- 
ity, as indeed it does. This however results in a local 
relaxation time which changes by many orders of mag- 
nitude over just a few particle spacings. This time scale 
gradient would be spread out via the dynamics of mode 
coupling theory, diffusing the excess mobility deeper into 
the bulk. In first approximation, MCT has been shown 
to correlate dynamics on the length scale £mct[I], the 
predicted divergence of which at T c will be broken by the 
emergence of activated events [25, 26J. A simple equation 
describing this gradient of relaxation rate incorporating 
activated events in the framework of Bhattacharyya et 
al.[2"5] emerges 



e M CT^ 2 r-\z) - (r-\z) - r^ cal {z)) = 0. 
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If we treat the mobility profile from purely activated 
dynamics, tZz „i(z), as a set of boundary conditions at 
z = and z = oo (equivalent to coarse graining on the 
length scale £,cr.r), then the relaxation time, decaying 
smoothly from r sur f oc e~ A °"' r f' Sc on the free surface to 
Tbuik oc e~ Abulk ' Sc <C r sur f in the bulk on a length scale 
£mct, would follow 
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The distance the excess mobility penetrates, z* , is found 
by comparing the magnitude of the two terms. 
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that Abuik/s c (T g ) w 40 at the glass transition. This gives 
a length scale z* Rj 20£mct which can be much larger 
than the bare dynamical correlation length. 

At a depth z < z* from the surface the 
particle motion would freeze out when r _1 (z) Rj 

p -A surf /s c (T l g °° al ) p -z/tMCT 



T bulk(Tg), 



or when 
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We find that for the thinnest films local glass transition 
temperature will be T™" = T° ur f as expected. T l g ocal 
will grow with distance from the free surface and be in- 
distinguishable from the bulk at a depth z* . 

The surface properties of glasses are important in many 
technological and biological contexts. Enhanced surface 
mobility is relevant for adhesion, friction, coatings, and 
nano-scale fabrication such as etching and lithography. 
Supercooled water at the surface of proteins acts to en- 
slave many protein motions[2T|l28|. Despite the influence 
of static surface perturbations, we feel that the ideal- 
ized treatment of the surface mobility of glasses presented 
here can help in understanding these phenomena. The 
enhanced mobility at free surfaces will allow phase trans- 
formations to occur at the surface that are kinetically 
impossible in bulk. De-vitrification often occurs at free 
glass surfaces, a fact of some importance in geology [29] 
and archaeology [30]. This observation is naturally ex- 
plained by the RFOT theory. 
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